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Abstract 

We present a numerical method of analyzing possibly singular incompressible 3D 
Euler flows using massively parallel high-resolution adaptively refined numerical 
simulations up to 8192 3 mesh points. Geometrical properties of Lagrangian 
vortex line segments are used in combination with analytical non-blowup criteria 
by Deng et al [Commun. PDE 31 (2006)] to reliably distinguish between singular 
and near-singular flow evolution. We then apply the presented technique to a 
class of high-symmetry initial conditions and present numerical evidence against 
the formation of a finite-time singularity in this case. 

Keywords: Euler equation, Existence, uniqueness and regularity theory, 
Vortex line geometry, Adaptive mesh refinement 



1. Introduction 

For now more than two centuries, the incompressible Navier-Stokes equa- 
tions have withstood the minds of mathematicians and physicists alike: The 
derivation of the nature of turbulence from the equations, as well as the global 
existence of smooth solutions is not known to date. The huge mathematical 
difficulties concerning the latter problem were recognized by its elevation to the 
status of "Millennium Prize Problem" by the Clay Mathematics Institute (see 
the official problem description by Fefferman [T], or review articles e.g. 
A proof of existence of global regular solutions to the Navier-Stokes equation is 
believed to entail the development of completely new methods for the analysis 
of partial differential equations. The absence of mathematical certainty for the 
Navier-Stokes equations may seem to leave the physicist in a somewhat embar- 
rassing position: The equation is known and well tested in application, but the 
existence of solutions is unclear in relevant cases. 

Yet, the actual impact of a supposed breakdown of solutions for the Navier- 
Stokes equations on physics of fluids is smaller than one might expect and 
appears like a mere technical detail on second thought. Singularities in the 
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Navier-Stokes equation would, if existent, appear on very small scales. Obvi- 
ously, continuum mechanics do not hold on these smallest scales and the break- 
down of the model equation would appear in a regime in which the model does 
not describe reality at any rate. Furthermore, the nature of supposed singu- 
larities for the Navier-Stokes equation is proven to be unphysical in nature, as 
it requires the existence of infinite momentum. Without external forcing, from 
smooth initial conditions and in the presence of friction, the occurrence of infi- 
nite momentum is impossible to justify physically. The impact of singularities 
is additionally limited by the fact that the space-time dimension of the singular 
region is proven to be less than or equal to one for the Navier-Stokes equations 

a- 

In the inviscid limit, the situation is quite the opposite. The incompressible 
Eulcr equations for ideal fluids, 

- U +u-Vu + Vp = 0, V-u = 0. (1) 
at 

appear to be of little physical significance in most applications, since friction is 
the dominating process on small scales. The ignorance regarding existence of 
global solutions is even larger for the inviscid case: The notion of weak solu- 
tions, which are well established for the Navier-Stokes equations since Leray jS], 
is unknown for the three-dimensional Euler equations. Nevertheless, the forma- 
tion of finite-time Euler singularities significantly concerns our understanding of 
(viscid) fluid dynamics. Euler singularities, if existent, would coincide with the 
development of large gradients in the velocity field. Since no friction is limiting 
the increase in velocity gradients, infinite momentum is not mandatory for a 
blowup of the Euler equations. The inviscid limit is, therefore, not only a mere 
description of ideal fluids, but explores the possibility of inherent dynamical 
processes beyond friction that limit the transition to smaller and smaller scales. 
This has immediate implications on the existence of a cut-off velocity in high 
Reynolds-number Navier-Stokes flows, leading to the slightly exaggerated ques- 
tion, quoting Constantin [6 : "Do we need Schrddingers equations to calculate 
the flow around a moving car? Or to predict tomorrow's weather?" For that 
reason, the problem of singularities for the Euler equations is of far greater im- 
portance to the physical understanding of fluids than the analogous problem for 
the Navier-Stokes equations. 

A similar argument is valid for turbulence. Today's phenomenological de- 
scription of turbulence (e.g. [3 [5]), which is built on the basis of the cele- 
brated theory by Kolmogorov [5HTT]. contains as a central point that, in the 
limit of vanishing viscosity, energy dissipation has to stay finite. This behavior 
could be explained by the formation of finite-time Euler singularities, as implied 
by Onsager's conjecture [12 . For three-dimensional incompressible flow, non- 
conservation of energy might be caused not only by viscosity but by missing 
regularity in the velocity field. Energy dissipation might occur, if the Holder 
continuity exponent is smaller than 1/3 for the velocity field. This conjecture 
was proven in terms of Besov spaces [T3l Ej • As a consequence, a mathemat- 
ical description of turbulence might be possible in terms of weak solutions for 
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the Euler equations, if smooth solutions gain enough roughness in finite time. 
Therefore, insight into the formation of finite-time singularities for the Euler 
equations could uncover a mechanism essential for the understanding of viscous 
turbulence. 

The search for finite-time singularities of the Euler equations has resulted 
in extensive literature, with many analytical results being relatively young. Es- 
pecially the advent of scientific computing has given research a new direction: 
Reports of numerical evidence supporting or denying the existence of finite-time 
singularities for the Euler equations are numerous (see e.g. [15] for a compiled 
list). 

As a now classical result, the blowup criterion of Beale et. al [TB] (BKM) 
connects the existence of solutions for the incompressible Euler equations in 
three dimensions to the critical accumulation of vorticity. More recently, geo- 
metric analysis of the flow [17l [18] has helped increasing insight into the process 
of vorticity growth. Among these geometric blowup criteria, theorems devel- 
oped by Deng et. al [H] may be seen as the first to be suitable for verification 
by direct numerical simulations. An approach along this way will be presented 
in this paper. 

Given the results of analytical considerations and the experience gained from 
numerical simulation of the Euler equations, certain scenarios are known to 
be possibly compatible with the analytic requirements of a finite-time blowup, 
namely the global notion of self-similar collapse to a point and the local pro- 
cess of vorticity accumulation by vorticity-strain coupling. It has been tried 
in the past to construct explicit initial conditions exploiting these scenarios to 
obtain numerical evidence for or against a finite-time singularity, with surpris- 
ingly inconsistent results. The major reason for this ambiguity is the critical 
dependence on extrapolation, which renders the identification of singular ver- 
sus near-singular behavior next to impossible by numerical means. The hopes 
are high that the situation is less vague when considering geometric properties 
of the flow, as mentioned above. We will present the application of such ge- 
ometric criteria to numerical data to sharpen the distinction between singular 
and near-singular flow evolution and identify the processes connected with this 
behavior. 

This paper is organized as follows: we first review the relevant geometric 
blowup-criteria that form the basis of our numeric method. This includes a 
motivation of how geometric properties such as curvature or spreading of La- 
grangian vortex line segments are connected to the accumulation of vorticity. 
We then describe the special class of high symmetry initial conditions used for 
our numerical experiments, discuss the implications of flow symmetries on the 
process of vorticity-strain coupling and introduce different vorticity profiles for 
vortex dodecapole initial conditions. Using this setup, details of the its imple- 
mentation for our massively parallel simulations with up to 8192 3 mesh points 
are given. Results are presented concerning the growth of vorticity and strain, 
the BKM-criterion and the geometric criteria. These findings act as numerical 
evidence against the formation of a finite-time singularity for this class of initial 
conditions. A conclusion and outlook summarize the paper. 
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2. Geometric blowup criteria 

Classical criteria for the development of a finite-time Euler singularity have in 
common that they focus on global features (such as certain norms of the velocity 
or the vorticity fields) or on point-wise Eulerian features (such as f2(t)) of the 
flow. This comes at the disadvantage of neglecting the structures and physical 
mechanisms of the flow evolution. A strategy to overcome such shortcomings 
was established by focusing more on geometrical properties and flow structures, 
such as vortex tubes or vortex lines. Starting with the works of Constantin et al. 
[T7JI20] and Cordoba and Fefferman [TS], some of these "geometric" criteria (e.g. 
p~9| EU [22]) have reached a phase where they allow direct verification of their 
assumptions with the help of numerical simulations. Special focus is placed on 
the criteria presented by Deng, Hou and Yu [IHIOSI], as the assumptions are in 
close reach for numerical simulations. They examine the Lagrangian evolution 
of vortex line segments and formulate a combined bound on velocity blowup 
and vortex segment collapse. 

One of the trivial consequences of the BKM theorem is the fact that no 
blowup can occur for the two dimensional Euler equations. Since the vorticity 
w(x, t) is bounded by the initial conditions ||u>o||.£°° for all times, a critical 
accumulation is impossible. This is a direct consequence of the vorticity pointing 
out of the plane of motion, therefore having the vortex-stretching term u> ■ 
Vu vanish everywhere. This may be interpreted as a motivation to focus on 
the behavior of the direction of vorticity, £ = u>/|u>| in the three-dimensional 
case. For 2D, £ is a constant in space and time (neglecting sign). In 3D, the 
consequences of the regularity of £ on the growth-rate of vorticity and ultimately 
of the applicability of BKM can be precisely stated. 

For the Euler equations, this was introduced by Constantin et al. |17j . They 
state, roughly, that for a smoothly directed vorticity in an 0(l)-region there may 
be no blowup in finite time as long as the velocity remains finite in this region. 
Even though this criterion takes into account the local structure of the flow and 
follows the evolution of vortex lines, the (global) bound on the velocity makes 
this theorem hard to apply in practice. Numerical simulations of the Euler 
equations give no evidence for the velocity to be uniformly bounded in time. 
This restriction on the velocity field is weakened in a similar criterion by Cordoba 
and Fefferman |18j . They consider vortex tubes with some properties concerning 
their regularity and a surrounding 0(1) region Q of the flow. From this it is 
possible to deduce, with the help of a milder assumption on the surrounding 
velocity, that the vortex tube cannot reach zero thickness in finite time. Even 
though the velocity field is no longer required to be uniformly bounded in time, 
the notion of "regular tube" of O(l) length is too restricting, compared to the 
experiences of numerical simulations. 

2.1. Regularity of vorticity direction along a vortex line 

Vortex lines of the three-dimensional incompressible Euler equations, defined 
as integral curves of the vorticity direction field, are transported with the flow. 
As a consequence, two points x and y on the same vortex line c(s) stay on 
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Figure 1: Two ways to apply theorem 1. Left: Choose y(t) such that it is far outside of the crit- 
ical region of maximum vorticity and monitor the behavior of f*(*} (V ■ $,) (c(s),t) ds. Right 

ry(t) 



For the position x(t) of maximum vorticity, choose y(t) such that ffj (V ■ £) (c(s), t) ds 
For a point-wise singularity, x(t) and y(i) must collapse in finite time. 
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the same vortex line for all times. Furthermore, as a direct implication of the 
solenoidality of the vorticity vector field, the vorticity flux through a vortex 
tube is the same for each cross-section. 

These two arguments may be combined to get a differential notion of the 
connection between the vorticity at two different points on the same vortex line. 
A simple consequence of the solenoidality of uj results in 

(€-v)M = -M(v-€). (2) 

Since for a vortex line c(s) it holds by definition that c(s) = £(c(s)), we have 
£• V = d/ds, where d/ds is the partial derivative in direction of the vortex line. 
Thus, integrating eq. |2| along the vortex line yields 



,(y(t),t)\ = \u(x(t),t)\ exp (- ^ V • £ds j . (3) 



Paraphrased, this means: The vorticity at two different points on the same 
vortex line is connected by the amount of converging or diverging of neighboring 
vortex lines along their interconnecting path. The more violent vortex lines 
converge around a vortex line, the faster the vorticity grows along that line. 

This finding was connected with BKM by Deng et al. [TH] to formulate a 
geometric blowup criterion. It is obvious from equation (|3| that the maximum 
vorticity f2(t) at a given time t can be estimated by the vorticity on its vortex 
line, as long as V • £ remains finite. In detail this means: 

Deng-Hou-Yu theorem 1: Let x(t) be a family of points such that for 
some Co > it holds \u)(x(t),t)\ > cnO(t). Assume that for all t G [0,T) there 
is another point y(t) on the same vortex line as x(t), such that the direction of 
vorticity £(x, t) = u>(x,t)/\cj(x,t)\ along the vortex line c(s) between x(i) and 
y(t) is well-defined. If we further assume that 



y(t) 

(V-€) (c(s),t) ds 

x(t) 



<C (4) 
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for some absolute constant C, and 



[ |w(y(i),t)|dt<oo J 



(5) 
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then there will be no blowup up to time T. 

It is immediately clear how this criterion can be applied to numerical sim- 
ulations: If the maximum vorticity £l(t) exhibits fast growth in time for which 
it is hard to decide whether it is a finite-time blowup compatible with BKM, 
instead one could monitor the vorticity outside the critical region, but on the 
same vortex line. If it remains well bounded, and V • £ along the vortex line 
does not scale critically in time, it is safe to deduce a non-critical growth of 
fl(t). This approach is sketched in Fig. [I] (left). 

However, due to the freedom of the choice of y(t) on the critical vortex line, 
theorem 1 may be employed in a different way to distinguish different scenarios 
for a finite-time singularity, as depicted in Fig. [I] (right). Suppose that the 
maximum vorticity Q(t) grows in a way compatible with BKM. Now, choose 
x(i) to be the position of maximum vorticity and define y(t) via 



for some constant C independent of the time t, where s denotes the arc-length 
parameter of the curve from x(i) to y(t). In words, choose y(t) on the same 
vortex line as x(£) such that the accumulation of tightening of nearby vortex 
lines is the same for every instance in time. This provides us with the ability to 
clearly distinguish between to separate cases of supposed blowup: 

(i) For every constant C, y(t) approaches x(i) in finite time to collapse to a 
single point. This would constitute the desired behavior for a point-wise 
singularity in the origin. 

(ii) If for any constant C, x(t) and y(t) stay well separated in time and do 
not collapse to a point, the whole vortex-line from x(t) to y(t) has to blow 
up in order to maintain critical growth in x(t). This scenario, however 
unlikely, is not ruled out analytically. 

The insight provided by the theorem could successfully be used as evidence ex- 
cluding a point-wise singularity for the considered initial conditions as presented 
below. 

2.2. Vortex line stretching and vorticity accumulation 

Vortex stretching is recognized as the mechanism responsible for the accu- 
mulation of vorticity. Revisited from a geometric point of view, vortex lines are 
transported with the flow, yet twist and turn due to vortex stretching. Since 
in the absence of dissipation vortex lines are unable to reconnect, the topologi- 
cal properties of vortex lines are fixed. A complex flow will therefore entangle, 
stretch and twist vortex lines in a non-trivial way. 



/ V • |d.s = C 

y X ( t ) 



(6) 
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Figure 2: Characterizing vortex line geometry in terms of \(Lt). A slightly curved vortex 
line with approximately parallel neighboring vortex lines (left) exhibits small X(Lt). Vortex 
lines with tightening neighboring vortex lines (center) or vortex lines with high curvature, in 
comparison to their length (right) have high \(Lt). 



The geometric equivalent of the vortex stretching term is the increase in 
length for a Lagrangian vortex line. It is possible to quantify this stretching 
and establish a sound connection to the vorticity dynamics of the flow. This 
in turn can then be used to reformulate blowup criteria in terms of geometric 
constraints on Lagrangian vortex lines. This section is meant to give an overview 
over this procedure to understand the implications of the second theorem of |19] . 
Details regarding the statements below are given therein. 

Consider a vortex line segment Lq at time t — and its Lagrangian image 
L t = X(Lo,£). Let f3, s be the arc length parameters of L t at times and 
t. Then, a direct implication of the vorticity transport formula describes the 
evolution of the absolute vorticity at a Lagrangian fluid element 

\u,(X(a,t),t)\ = £(X(a,t),t) • V„X(a,t) • C (a)|« («)| (7) 

= 0l«o(a)|, (8) 

meaning that the local stretching of the length of a vortex line segment that 
is transported with the flow is equivalent to the growth of vorticity at the 
corresponding transported fluid element. 

This result can be transformed into a bound for the length of a vortex line 
by the vorticity along this line. Denote with l(t) the length of the vortex line 
segment L t at time t and define with 

n £ (i):=||«(.,i)IU~(2iO (9) 

the maximum vorticity on the vortex line segment. Furthermore, let 

M{t) := max(||V • t\\ L ~(L t ), INU~(£ t )) (10) 

be the quantity of vortex line convergence V • £ and vortex line curvature k, and 
define X(L t ) := M{t)l(t). Then, the relative increase of the length of the vortex 
line segment in a time interval, l(t)/l(0), is bounded as 

e -X(L t) Wl < m <e HL ) S m Ml) 

n,(o)-i(o)- e fij(o) • [Li> 
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Figure 3: Visualization of the evolution of a vortex line segment in theorem 2: At a late time 
t'z > t\ the segment Lt 2 has to be included in the Lagrangian evolution of of Lt 1 , but is free 
to be just a fraction of it. 



This result, which is a slight modification of a lemma in |19j . may be viewed 
in its own right: The relative increase in length along a time-interval is bounded 
by the vorticity increase and a factor exp(±A(L t )). Thus, X(L t ) is a dimension- 
less number, characterizing the geometric "tameness" of the vortex line filament. 

As depicted in Fig. [2j a vortex line segment has a huge X(L t ), if its maximum 
curvature is large, relative to its length (the segment is "kinked" instead of 
"curved" ) , or if the surrounding vortex lines collapse to the considered segment 
in at least a point (the surrounding is "tightening" instead of "parallel"). A 
relatively unbent vortex line segment with approximately parallel neighboring 
vortex lines possesses a low value of X(L t ). This quantifies the constricted notion 
of "relatively straight" and "smoothly directed" given in [IT] in a sharper way. 

2.3. Lagrangian evolution of vortex line segments 

Connecting the stretching process stated above with the Lagrangian accu- 
mulation of vorticity, 

^M = [(S€)-e]M, (12) 

where S = 1/2 (Vu + Vu T ) is the strain tensor, and noting that the curvature 
k of the vortex line L t fulfills 

— (13) 

with n = L t /\L t \ being the unit normal vector of the vortex line, the Lagrangian 
evolution of vortex line stretching, by inserting eq. (j8|), becomes 

D ( ds\ d . . . f ds\ 

= ^( u '0-«(u-n) — . (14) 



Dt \df3 J dp y ' \d/3 

At this point it becomes obvious how the process of vortex line stretching in- 
teracts with the velocity in two distinct ways: The velocity in direction of the 
vortex line elongates the segment by drawing it out, while a part of the velocity 
normal to the vortex line increases the segment's length by enlarging its curves. 
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Integrating (14) along the vortex line (from fti to P2) and over time (from 
to t) results in 

lit) < 1(0) + f [U s (t) + X(t)U 71 (t)} dr , (15) 
Jo 

for 

U ( (t):= max |(u • £)(x, t) - (u • £)(y, t)\ 
x,yeL t 

U n (t) := max |u • n| 

Instead of starting the above reasoning at time t = 0, the results are identical for 
a later time < t\ < t. This result may be understood as an upper bound for 
vortex line stretching in terms of velocity and vortex line geometry. In conjunc- 
tion with the connection between length increase and vorticity amplification, 



given in equation (11 1, one arrives at 

r 

1(0) 



fi,(t) < fi,(0)e A(it) 



1 + 77^ / (^(r) + Hr)U n (r))dr 



This is an inequality for the control of growth rate of the vorticity by geometric 
flow properties. From this estimate, by combining it with BKM to distinguish 
critical from sub-critical vorticity growth, the central non-blowup criterion of 
[19] is derived: 

Deng-Hou-Yu theorem 2: Assume there is a family of vortex line seg- 
ments L t and Tq £ [0, T), such that L t2 C X(L tl , t±, £2) for all Tq < t\ < t^ < T . 
We also assume that VL(t) is monotonically increasing and ||a;(i)||j,oo(£ t ) > 
coQ(t) for some cq > when t is sufficiently close to T. Furthermore, we 
assume that 

(i) A(L t ) < Co, 

(ii) l(t) >(T- t) B for some B G (0, 1). 

(iii) U^t) + U n (t)X(L t ) < (T - t)~ A for some A < 1 - B 

Then there will be no blowup in the 3D incompressible Euler flow up to time T. 



Here, a(t) < b(t) means there exists a constant c 6 K such that \a(t) \ < c \ b(t) \ 
(and accordingly for a(t) > b(t)). The choice of Lagrangian vortex segments 
is sketched in Fig. [3j 

It should be noted that theorem 2 again includes assumptions on the di- 
mensionless number \(L t ). Especially assumption ^ poses a uniform bound 
in time for X(L t ). This translates to words as the process of "zooming in" to 
the location of maximum vorticity in order to keep the considered vortex line 
segment relatively straight in comparison to its length. The assumed accompa- 
nying collapse in length to keep X(L t ) bounded is then linked in its growth rate 
to the blowup of the velocity components. 

It is worth mentioning that the above presented criterion, even though it 
is obviously inspired by the classical geometric criteria, still differs in crucial 
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aspects. The posed assumptions are purely local and restricted to the geometry 
of a single critical vortex line filament. Assumptions on the velocity do not, 
in contrast to Constantin et al. pT], impose a uniform bound (which is not 
observed in simulations), but allow for a finite-time blowup of velocity, strictly 
connected in its growth rate to the geometrical evolution of the filament. The 
vortex line segment itself is not assumed to be of O(l) length (as in [18 J or to 
be contained in an 0(l)-region (which, again, was not observed in simulations). 
These aspects in combination render it a promising theorem to be directly tested 
by numerical simulations. 

2.4. Scenarios for finite-time singularities 

It has been established that a singularity of the Euler equations in finite time 
necessitates rapid accumulation of vorticity. Locally, vorticity-strain coupling is 
identified as the mechanism for nonlinear amplification in finite time. Globally, 
the notion and possibility of self-similar or locally self-similar collapse to a point 
is introduced. These aspects serve as a basis for the construction of initial 
conditions suitable for the possible formation of a finite-time singularity. 

Due to the incompressibility condition, the trace of S vanishes and due to the 
symmetry of S its eigenvalues A; are real and the corresponding eigenvectors 
Vj are orthogonal. Thus Ai, the biggest eigenvalue, fulfills Ai > in regions 



with non- vanishing strain. As a direct consequence of equation (12 1, if Ai is 
proportional to the vorticity in direction Vj then the vorticity growth would be 
compatible with a finite-time singularity according to BKM. This coupling of 
the strain to the vorticity is crucial. If the strain rate is constant instead, the 
growth in vorticity is merely exponential. Several cases have been suggested in 
which this mechanism of coupling may take place. 

It should be noted that, contrary to expectations, the vorticity-strain cou- 
pling does not readily appear in nature. One would expect a tendency of the 
vorticity to align itself to the eigenvector of the strain tensor with the largest 
eigenvalue all by itself, since the parallel component is amplified, while the or- 
thogonal components are damped or stay nearly constant. Nevertheless, for 
viscid turbulent flows, quite a different behavior is observed both in numerical 
simulations and experiments: The vorticity is most likely to align to the in- 
termediate eigenvector of the strain tensor 24 26 and similarly for the Euler 
equation |27j . One can therefore expect that functional vorticity-strain coupling 
is inherently unstable. The process has to be designed "artificially" by choosing 
suitable initial conditions. 

A more precise notion of the process of vorticity alignment in turbulent 
flows is given by Hamlington, Schumacher and Dahm [28]. They distinguish, 
evaluating the Biot-Savart law numerically, between strain induced locally by the 
immediate neighborhood and globally by long-range interaction. For turbulent 
flows, they observe a most likely alignment of the vorticity to the most positive 
eigenvector of the global strain. Taking into account also the local strain restores 
the alignment of vorticity to the intermediate eigenvector. Contrary to this, for 
the successful emergence of a finite-time singularity of the Euler equations in a 
point-wise sense, vorticity-strain coupling should be induced by the local strain 
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Figure 4: Effects of planes of refiectional symmetry on intersecting vortex tubes: In the 
symmetry plane, the vorticity u> is normal and the strain tensor possesses a parallel eigenvector 
with corresponding eigenvalue 5 ZZ . The curvature n in the symmetry plane has to increase 
in order to support critical vorticity-strain coupling. 

to support a collapse to a point. A numerical application of this technique from 
|28j to Euler blowup simulations will be presented elsewhere. 

3. Initial conditions 

Different initial conditions were introduced and subsequently improved or 
refined to construct flows with prolonged intervals of vorticity-strain coupling, 
and it seems natural to search for techniques to artificially keep the coupling 
existent. One such technique is the introduction of symmetries to the flow. 
Early examples such as the Taylor-Green vortex or Kerr's initial conditions 
|30H32j are employing such symmetries. 

3.1. Refiectional symmetries 

Consider the plane z = to be a plane of refiectional symmetry, as shown 
in Fig. [4j defined by 



for the velocity vector field, which leads to u z = in the plane of symmetry. 
Accordingly, the vorticity obeys 



u x (x,y,z) 
u y (x,y,z) 
u z (x,y,z) 



u x {x,y,-z) 
Uy(x,y,-z) 
-u z (x,y, -: 



u x [x,y,z) = -u x {x,y,-. 
uj y (x,y,z) = -u v (x,y,-. 
uj z (x,y,z) =u z (x,y,-z) 



z) 
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and in particular lu x = oj v = or oj = u z e z in the plane of symmetry. Due to 
these properties, the strain tensor has S xz — S yz = S zx — S zy = 0, or 



S=\S Xy Syy . (16) 




It immediately follows that the eigenvector corresponding to the eigenvalue S zz 
is directed normally to the symmetry plane, and the vorticity vector is aligned 
to it. Note that this is the sole consequence of the reflectional symmetry and is 
in no way influenced by the flow. A vortex tube normal to the symmetry plane 
therefore seems like a natural candidate for critical accumulation of vorticity 
by means of vorticity-strain coupling: All that is needed is a sufficiently long 
period of time in which S zz ~ uj z at one point of the symmetry plane. 

This possibility has been analyzed by Pelz [53]: Taking into account only 
the zz-component, the strain tensor in the plane of symmetry is given by the 
Biot-Savart law as 

Szz = i I {{x ~ x>y{x ' ] ~ {y ~ feS^' ■ (17) 



Equation (17) shows that S zz in the plane of symmetry does not scale with u z , 
but does instead depend on uj x and u> y , which are both equal to zero in the z — 
plane. Yet, in close proximity to the plane, ui x and uj y may grow, depending 
on the curvature of the vortex line intersecting the symmetry plane: A Taylor 
expansion around the z — plane yields: 



duii 

h IT" 
oz 



= hKiUj z (z = 0) 



=o 



for small h up to first order, for Ki — ftrij. Therefore, if the curvature is huge 
close to the plane of symmetry, with i g {x, y} approximately scales like 
u> z and thus S zz may scale with uj z too. However, for this to happen, we need 
k ~ 1/h. As a matter of fact, the dimensionless number Kih plays a similar 



role as the characteristic geometric number X(t) introduced in section 2.2 For 
S zz to blow up like lu z , the curvature has to increase in a way to counter the 
shrinking of the length scale h. 

On the other hand, the axial strain S zz stretches the vortex tube in z- 
direction. This counteracts any increase in curvature to a certain degree. More 
precisely, the Lagrangian evolution of the curvature components k x and K y were 
calculated in [55] : 

^y^ X = i^XX ~ 2Szz) K X + SyyKy + 9 Z S yz 

~jyL^V — ^xx^x i^yy zz )K y -\- d z S xz . 

The axial strain S zz diminishes both n x and n y . 
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These counteracting processes of vortex line geometry are by no means ana- 
lytically exact, since all long-range interactions have been ignored. Nevertheless, 
they constitute an intrinsic resistance of a single vortex line to "self-stretch" in 
a critical way. The argument may be readily translated to the case of (per- 
turbed) anti-parallel vortex tubes: Since no other components of vorticity are 
introduced, S zz still only depends on ui x and lu v which in turn rely on high 
curvature to scale like lo z close to the plane of symmetry. 

One way to counter this is to induce the axial strain by neighboring tubes 
instead of relying on a sufficiently large kink. This will be presented in the 
following section by introducing additional rotational symmetry. 

3.2. High symmetry initial conditions 

One notable high-symmetry flow was introduced by Kida [34 and has sub- 
sequently been used extensively to probe a possible Euler blowup numerically 
[3"5H3"7] or analytically (e.g. [55] ) as well as study the onset of turbulence (e.g. 
[39]). The Kida-Pelz flow has a three- fold rotational symmetry about the diag- 
onal and a reflectional symmetry about all three Cartesian planes. Flows with 
these two properties are termed as invariant under the full octahedral group [33] . 
The Euler (and Navier-Stokes) equations preserve the Kida-Pelz symmetries. In 
general one can write these initial conditions as 

v{x, y, z) — aimn s\n{lx) cos(my) cos(nz) (18) 

Z,m,n 

u = (u x ,u y , u z ) T = (v(x, y, z),v(y, z, x),v(z, x, y)) . (19) 

This means, for a computational domain spanning the interval [0, n] in all three 
dimensions, that the normal component of the velocity field is anti-symmetric 
under reflection at the Cartesian planes while the tangential components are 
symmetric. Combining this with the three-fold rotational symmetry adds up to 
a total memory saving factor of 1 /24. 

On the same time there is reason to hope that these rather artificial sym- 
metries encourage singular behavior if the initial conditions are constructed 
accordingly. When assuming a localized vortex tube intersecting the symmetry 
plane normally, as depicted in Fig. [5] its mirror images result in a total of 
six pairs of anti-parallel vortex tubes. It has been proposed by Pelz [33] that 
the strain induced by the rotational images of each tube, assuming a velocity 
field supporting a collapse to the origin, may lead to the desired vorticity-strain 
coupling without being subject to the counteraction of strain and curvature in 
the planes of symmetry. Provided that the vortex dodecapole retains its shape 
during collapse, this scenario could lead to a finite-time collapse to the origin. 

3.2.1. Vortex dodecapole initial conditions 

One form of these initial conditions is based on the idea to already start with 
six dipoles consisting of vortices of a designated vorticity profile. An example 
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Figure 5: Left: Vorticity profile of one tube of the 12-tube initial condition (with A = 20). 
Right: Volume plot of the vorticity for the whole domain. 



is the vortex dodecapole initial condition [40] with a vorticity profile given by 



A 



1 — exp 



-e 2 log(2)— exp 
3r 



- 1 



for r < 



for r > - 
- 3 



(20) 



where r denotes the distance to the tube's center line. The vorticity decreases 
with increasing distance to the center line and is strictly zero for r > 2/3. 
Thus, the vortex has compact support in the r-tp-pl&ne while still being smooth. 
Fig. [5] (left) displays the vorticity profile given above, Fig. [5] (right) shows 
the whole dodecapole. Only one octant, i.e. three vortices, are simulated due 
to symmetry. For the reasons lined out in the previous sections, this kind 
of dodecapole appears to very promising in terms of developing a finite-time 
singularity in the origin because of reciprocal strain of the mirror tubes. It 
is furthermore susceptible to the analysis by the presented geometric blowup 
criteria. Most of the diagnostics in this paper are therefore performed on flows 
which are based on initial conditions of this type. 



3.2.2. Lamb- dodecapole initial conditions 

For the Euler equations, a single stretch-free axisymmetric vortex may have 
arbitrary radial dependence for the vorticity to remain stationary in time. The 
same is not true for vortex dipoles: An isolated vortex dipole propagating 
through the domain does not preserve its shape. This so-called vortex shed- 
ding is believed to influence and possibly suppress a self- amplifying behavior 

EH- 

There are exact form-preserving dipole solutions of the 2-dimensional Euler 
equations which may be used to construct initial conditions that do not exhibit 
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Figure 6: Comparison of vorticity profile for a Lamb dipole to the simple dipole. Left: 
Lamb dipole used in the Lamb-dodecapole initial conditions. Right: Dipole used in vortex 
dodecapole initial conditions. Both are scaled to fit in amplitude and size. 



vortex shedding. The most famous is the Lamb-dipole introduced by Lamb 
Following [43] it is defined by 



■lUk^MA sin(0) for r < a 



J (ka) 



(21) 







for r > 



where k is chosen such that ka is the first zero of J\, i.e. ka « 3.8317. In Fig. [6] 
this vorticity distribution is compared to the vorticity profile given in equation 
(20). Note that even though the distribution of vorticity for the Lamb dipole 



appears to be less sharp than for the simple profile, it is not differentiable at 
r = a, as can be seen in equation (21). Because of this, strictly speaking, the 



Lamb dipole is an improper candidate for the search for finite-time singularities. 
This issue is usually overcome by smoothing high frequency components in order 
to smear out the discontinuity in the gradient of the vorticity. 

Orlandi and Carnevale [H] where the first to use Lamb dipoles to construct 
a colliding pair of dipoles, observing a rapid amplification of vorticity for a 
period of time, with a slowing growth at later times due to either depletion of 
nonlinearity [44] or lack of resolution. The Lamb dipole is used in a similar 
manner in the context of this paper to form a Lamb dodecapole analogous to 
vortex dodecapole initial conditions presented above. 



4. Numerical experiment 

Along the lines of the above presented mechanism for a finite-time singular- 
ity for the Euler equations, a number of numerical simulations were performed 
in the last decade to act as evidence for or against a blowup. Beginning in the 
early 80s of the last century, numerous numerical simulations with a variety of 
methods and schemes and differing initial conditions where tested, from Pade 
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approximants [3SJ , vortex-segment methods [301 |37J and vortex filament mod- 
els [31], to projection methods 03], pseudo spectral simulations [3"21 1571 ISO! 1ST] . 
Chebychew codes [30] and adaptive mesh refinement j27j[S2]. Despite ever grow- 
ing resolution, from 128 3 [30] up to 4096 3 mesh-points of adaptive simulations 
[30] or latest pseudo-spectral codes [S3], results are often inconclusive or even 
conflicting. 

In this paper The need for high resolution is met via massively parallel 
adaptive mesh refinement for the results presented. This high resolution data is 
then analyzed on the basis of the geometric blowup criteria presented in section 
® 

4--1- Computational Framework 

All numerical simulations throughout this work were conducted using the 
recently developed framework racoon III (refined adaptive computations with 
object-oriented numerics, based on [S3])- Its key feature is the integration of 
partial differential equations on adaptive grids on massively parallel distributed 
computers. Most Euler blowup scenarios feature extremely localized structures 
with steep gradients, where a fixed mesh would under-resolve the crucial parts 
while wasting resources on the less important ones. 

Adaptive mesh refinement increases the locally available resolution, but 
comes at the cost of additional computational overhead. It complicates the 
framework in several ways. Most importantly, it restricts the choice of numeri- 
cal schemes to comparatively simple low order finite difference or finite volume 
variants. A direct comparison to high accuracy pseudo-spectral simulation was 
made in [30] for the case of Euler equations. It was found that a resolution 
approximately 1.3 times higher is needed to reach a comparable accuracy for 
the adaptively refined code. 

To decide which regions are to be refined, the norm of the gradient of velocity, 
||Vu(x, t)\\ is compared to a threshold. If the block is flagged as being under- 
resolved, it is bisected into 2 d child blocks that are redistributed among the 
available nodes. The resolution of the parent block is thus effectively doubled. 
The opposite happens for blocks that are over-resolved: 2 d blocks are merged 
into one, the resolution at this location is halved. With this procedure, the 
grid is constantly changing and adapting to the simulation, as shown in Fig. [7j 
allowing high resolution at critical locations but not wasting any resources for 
the rest. 

Since communication between different nodes is the smallest bottleneck due 
to limited bandwidth and high latency, it is advantageous to arrange the blocks 
in a way that physically close blocks are located on the same node. Even if this 
seems to be pretty straight-forward for normal grids, it poses a larger problem 
for adaptive grids with different resolutions. In racoon III, blocks are distributed 
along a space-filling Hilbert curve, as sketched in Fig. [S] 

This ensures that proximate blocks are located on the same node even if 
the grid is not fixed. Currently racoon III uses a slightly different approach, 
using independent Hilbert curves for each level, since inter-level communication 
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Figure 7: Refinement criterion for the simulation of the Euler equations with racoon III. 
Regions with a large value for ||Vu|| are resolved higher. Each square represents a block with 
16 3 cells. Shown is the absolute vorticity for a cross-section of one vortex tube. 
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Figure 8: Adaptive mesh refinement and dynamic load balancing. The workload is distributed 
among different processors along a space-filling Hilbert curve. 
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Figure 9: Mixed weak and hard scaling for the framework racoon for a hyperbolic test problem. 
The scaling is close to linear for up to 262144 cores. 



is the most frequent type of communication for common problems. Every time 
the grid changes when adapting to the current situation, the Hilbert curve is 
recalculated, as is the workload for each node. If an imbalance is detected, the 
blocks are redistributed along the curve, each node getting roughly the same 
amount of blocks. 

The numerical scheme consists of a strong stability preserving third order 
Runge-Kutta |55| time integrator combined with a third order shock-capturing 
CWENO scheme [56] to reduce oscillations in the presence of strong gradients. 
The integrated equation is the vorticity formulation of the Euler equations, 

d 

— u + V x (V(u® u)) = , (22) 

employing a vector potential formulation AA = — u) with u = V x A to ensure 
solenoidality of the vorticity vector field ui. The associated Poisson equation 
is solved with a second order parallel and adaptive multigrid algorithm via a 
full approximation scheme to account for the non-uniform meshes, with flux 
correction at refinement interfaces. Interpolation on the coarse-fine interfaces 
is done in ui for the whole scheme to ensure the highest possible accuracy in 
the critical variable. Passive tracer particles are injected into the flow for the 
tracking of Lagrangian vortex line segments. The above third order Runge- 
Kutta is also used for the time integration of the tracer particles and the space 
integration of vortex lines. 

The overall scaling for racoon is depicted in Fig. [9] for a hyperbolic test 
prolem (compressible MHD. It was measured on the BlucGcnc/P machine at 
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Forschungszentrum Jiilich with a total number of 294912 cores. For a combina- 
tion of weak and hard scaling, the performance is close to linear up to 262144 
cores, the maximum number tested. The elliptical problems encountered when 
simulating the Euler equations (velocity projection or calculation of the vector 
potential) are more difficult to parallelize than the hyperbolic advection term 
due to their inherent non-local nature. For each timestep, information trav- 
els only fractions of the grid spacing in the advection step, but through the 
whole domain when enforcing the incompressibility. This behavior is necessar- 
ily reflected by the demands on communication between processes in massively 
parallel simulations. With inclusion of the multigrid algorithm, the scaling is 
efficient only up to 131072 cores. 

4-2. Evolution of the flow 

This section is devoted to the visible results of the actual simulation of the 
presented vortex dodecapole configurations. For this purpose, the CWENO vec- 
tor potential formulation is used in conjunction with adaptively refined meshes 
for simulations with a resolutions of up to 8192 3 effective grid points, taking 
into account the increase in resolution due to the high symmetry of the initial 
conditions. Both the Lamb-dodecapole and the vortex dodecapole are used as 
initial conditions. 



4-2.1. Vortex dodecapole 

The vortex dodecapole was chosen as a prototype for the class of vortex 
dodecapole initial conditions. Its main features are a smooth vorticity profile 
with compact support and straight, unperturbed initial vortex tubes. 



Pictured in Fig. 10 is the evolution in time for the vortex dodecapole initial 



conditions. Shown are isosurfaces of the absolute vorticity |w(x, t)\ at 75% of the 
peak vorticity for different times. Due to the high symmetry, only one octant of 
the computational domain is simulated. The figures therefore depict only one 
half of a vortex tube, with twelve similar tubes in the total domain. The initial 
phase of the development is depicted in the first two sub-figures: The initially 
straight tube gets slightly stretched due to interaction with the neighboring 
tubes. In the third frame, the well-known flattening is in progress. The last 
three pictures present the final stage of the flow, where the tip of the sheet rolls 
up and forms a secondary vortex sheet. In the final figure, the secondary sheet 
exceeds the original sheet in length. Its tip gets drawn out of the collapsing 
region. 

The appearance of the roll-up and the secondary vortex sheet are a first 
evidence against a locally self-similar amplification and collapse to a point: The 
initially round vortex tubes are severely deformed and do not resemble their 
initial configuration in shape. Furthermore, the possibility that the formation 
of a roll-up may lead to the emergence of a tube-like structure which again form 
a dodecapole arrangement is clearly conflicting the numerical evidence. 
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Figure 10: Evolution of the vortex dodecapole. Pictured are isosurfaccs of the absolute 
vorticity, |o)(3:,t)| at 75% of the peak vorticity. Only one of twelve tubes is shown. The 
flattening of the vortex tube is followed by a roll-up. The developing secondary sheet finally 
exceeds the original sheet in size. All pictures are from run amrl. 
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Figure 11: Evolution of the Lamb dodecapole. Pictured are isosurfaces of the absolute vor- 
ticity, \uj(x, t)\ at 75% of the peak vorticity. Again, only one of twelve tubes is shown. As 
before, the vortex tube is flattens, followed by a roll-up. A secondary vortex sheet develops 
and gets drawn out of the center region. All pictures are from run lamb. 
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4-2.2. Lamb dodecapole 

The Lamb dodecapole initial conditions are motivated by the fact that each 
Lamb dipole in itself is an exact and invariant solution to the Euler equations. It 
was therefore anticipated in |41j that a more complex setup consisting of Lamb 
dipoles will exhibit considerably less core deformation for the vortex tubes. If 
this assumption would be met, the dodecapole arrangement could lead to the 
formation of a locally self-similar blowup scenario: The Lamb-dipoles would 
approach and amplify each other, but, without core deformation, stay in their 
relative alignment and shape. The ever-decreasing length-scale would result in 
a point-wise collapse to the origin. 

As shown in Fig. this scenario is not observed in the numerical simu- 
lation. The initial tubes are deformed severely in the course of the simulation. 
Vortex core deformation is not prevented. This is hardly surprising, since the 
vortex dodecapole relies on strain imposed by the rotational images of the tube 
by design, while the Lamb dipole configuration only prevents deformation by the 
reflectional image. Due to the initially close proximity of all twelve vortex tubes 
and the short timescale of the evolution, deformation induced by the reflectional 
partner seems to be negligible, regardless of the actual vorticity profile of the 
tubes. 

Altogether, the evolution of the vortex tubes for the Lamb case resembles 
the above presented vortex dodecapole flow: An initial flattening of the tubes 
is followed by a roll-up. The emerging secondary vortex sheet gets drawn out 
and finally exceeds the original sheet in length. Due to the overall similarity 
of both flows it seems safe to deduce that the topological flow evolution only 
weakly depends on the precise vorticity profile. This may be seen as motivation 
to transfer the results for just one particular initial condition to the whole class 
of vortex dodecapole flows. 

4-2.3. Comparison and conclusion 

The results of the previous two sections lead to the conclusion that no re- 
markable differences exist in the overall properties of the flow. In Fig. [l2j a 
direct comparison between low resolution runs (1024 3 ) for the simple and the 
Lamb vorticity profile are shown for a late time to reveal the details of the differ- 
ences. Most of the large-scale structures are identical for both flows. The initial 
shape of the Lamb profile is responsible for the formation of a less sharp roll-up 
of the vortex sheet and the accumulation of secondary vorticity inside the kink. 
Furthermore, the trailing vortex sheet, which is an artifact of the collapse of the 
vortex dipoles to the center, is considerably stronger for the Lamb dipoles. 

Since, additionally, the core deformation is not effectively prevented in the 
Lamb case, these arguments were the reason that all high resolution runs and all 
geometric diagnostics were performed for the simple, smooth dodecapole initial 
conditions. 

4-3. Accumulation of vorticity and strain 

The vortex dodecapole is designed to be a violent initial condition with rapid 
accumulation of vorticity. Unlike e.g. the Taylor-Green vortex or Kerr's initial 
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Figure 12: Direct comparison between the Lamb and the simple vorticity profile at late time. 
The vorticity in a slice near the plane of symmetry, z = 0.1, is pictured. The Lamb dodecapole 
(top) exhibits a more pronounced trailing vortex sheet near the symmetry plane. This effect 
is considerably smaller for the simple vorticity profile (bottom). 



conditions, no sustained phase of flow evolution has to be awaited for the critical 
structures to form. Thus, vorticity accumulation sets in immediately. 

The overall vorticity amplification from initially £7(0) = 20 to finally Q(i) > 
10 5 exceeds a factor of 500. The location of the maximum vorticity follows the 
tip of the vortex sheet show above, and is located at the intersection of the 
vortex sheets when the roll-up begins to form. The growth of the maximum 
of the norm of the strain, behaves in a similar manner as the peak 

vorticity, which is about \\S(-, 0)||l°° ~ 12.4 initially and grows by more than 
two orders of magnitude in the course of the simulation. 

The BKM-criterion implies that the growth in time of Q(t) has to fulfill 
fi(i) « 1/(T — t) 1 with 7 > 1 to be compatible with a finite-time singularity. A 
plot of (by assuming j = 1) is pictured in Fig. 13 At small times i, this 

graph looks straight, but the growth rate changes at least twice in the evolution 
of the flow. This can be explained by competing maxima in \u}\ overtaking the 
original fi(t), thus changing the growth rate at different stages. Nevertheless, 
at no time the vorticity looks as though saturating, and in the latest stage of 
development suggests a blowup time of T sa 0.72. 

Numerical data of this kind has been interpreted as evidence in favor of 
the formation of a finite-time singularity before. Yet, even though the plot [13] 
is rather suggestive, the growth may as well be fitted to some fast (double) 
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Figure 13: Evolution of 1/Q(t) in time. This mode of plotting suggests a growth of fi(i) RJ 
1/(T - t)T with 7 = 1 and a blowup time T sa 0.72. 

exponential growth. 

If.. If.. Geometry of the critical vortex line 

It was stated by theorem 1 of |19j that a blowup of vorticity in any point x 
is impossible as long as for some y on the same vortex line uncritical growth of 
vorticity is observed and along the vortex line connecting x to y the integral of 
V • £ remains bounded. As lined out above, we re-interpreted this statement as: 
Supposing there is singular behavior of the maximum vorticity does the 

flow allow for a point-wise blowup or is there a blowup of a finite, non- vanishing 
vortex line segment? 

Numerically, this test was implemented as follows: 

• At each timestep, identify the point of maximum vorticity as x(t). 

• Follow the vorticity direction vector field while integrating V • £ along the 
path. This is done with a third-order Runge-Kutta integrator in space. 

• As soon as the integrated quantity exceeds the threshold C, identify the 
current location on the vortex line as y(t). To increase precision, the 
endpoint is found via bisection. 

• Geometric properties and diagnostics for the vortex line segment are cal- 
culated, especially its length and \u)(y(t),t)\ to distinguish the cases in- 
troduced above. 

This procedure is carried out for the whole time interval, as long as the simula- 
tion is well resolved. The constant C is chosen in a reasonable way to achieve 
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Figure 14: Left: Length of the vortex lines starting at position x of maximum vorticity for 
constant C = f v V • £ds. Right: Vorticity at the endpoint y of these vortex line. Once 
satiated, the growth rate is the same for all y. 



a length of the vortex line segment that fits into the computational domain in 
the beginning of the simulation, but is still well resolved at the chosen resolu- 
tion at later times. Hence, the whole vortex line segment is resolved reliably 
throughout the simulation. 

The results for the vortex dodecapole initial conditions are presented in Fig. 



14 for different constants C € {0.25, 0.5, 1, 2}. Initially, the vortex line segments 



do not accumulate enough V • so that the length is bounded by the size of the 
computational domain (x G [0,7r] 3 ). At some point, depending on the value of 
C, the threshold is reached and the length of the vortex line segment decreases. 
Yet, for all considered cases of C, the length does not collapse to a point, but 
saturates at early times without approaching l{t) = 0. This behavior appears 
to be stable up to the latest time of the simulation. The final length of the 
vortex line segments is at least 0.3 for the smallest case of C (C — 0.25), which 
is still well resolved with at least 200 Ax. This result, therefore, is a numerical 
evidence against a point-wise blowup for the vortex dodecapole class of initial 
conditions. This is in concordance with the estimate in |19j . 

Monitoring the development of u>(y(t), t) yields, as shown in Fig. 14 (right), 
a similar growth rate for the accumulation of vorticity at the endpoint as for 
the beginning of the vortex line segment. This is hardly surprising, since by 
construction a constant value for C directly links the growth rates of |u;(x(i), t)\ 
to \uj(y(t),t)\. Nevertheless, a numerical verification of this analytic equality 
may be seen as a confirmation that the observed growth rate of \u>(x(t),t)\ is 
by no means a numerical artifact in an isolated small area, but is reproduced 
at points far away from the critical region, which appear to be well-behaved 
at first view. The possibly critical growth in the perspective of BKM is, thus, 
confirmed by the global flow. 

Furthermore, since for a large portion of the simulation the distance l(t) 
is approximately constant, this could possibly be seen as an evidence for the 
existence of a non-vanishing vortex line segment that blows up in every point. 
The popular scenario of a collapse to a single point, on the other hand, is clearly 
conflicting the numerical evidence. The discovery of a possibly critical vortex 
line segment in the vortex dodecapole flow, however, is to be handled with 
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care, since distinguishing between critical and sub-critical blowup of the whole 
segment is in no way more conclusive than distinguishing between critical and 
sub-critical growth of f2(t). Thus, learning from the lesson taught by 25 years 
of numerically testing BKM, this should not be interpreted as clear evidence in 
favor of a finite-time singularity. 

4-5. Lagrangian evolution of the critical vortex line segments 

The geometric properties of Lagrangian vortex line segments, especially their 
curvature k and the tightening of their surroundings V • £ have been established 
as revealing parameters in understanding the nature of rapid accumulation of 
vorticity in Euler flows and a sound connection to singular behavior is made 
through theorem 2. The ambition here is to utilize these geometric properties, 
monitored in a numerical simulation, as more reliable means of distinguishing 
between a finite-time singularity and a mere fast accumulation of vorticity. 

Despite high hopes from an analytical point of view that these considera- 
tions will shed light on the nature of vorticity accumulation, numerical results 
observing geometrical properties of Lagrangian vortex filaments are scarce. This 
is primarily due to the fact that Eulerian quantities such as Q(t) are readily 
trackable in post-processing, while monitoring the Lagrangian evolution requires 
additional computational effort. On top of that, the geometry of integral curves 
at an instance in time, though in principle computable in post-processing, as 
well as derived quantities such as their convergence and curvature, are quite 
inaccessible in comparison to simple Eulerian criteria. 

This section is devoted to the presentation of results concerning the assump- 
tions of theorem 2 of [TH] for the vortex dodecapole initial conditions. Quite 
similar to the first theorem, there is considerable freedom in the choice of the 
involved quantities. The strategy we chose in the context of this paper is as 
follows: 

• Identify the Lagrangian fluid element a, which will contain the maximum 
of vorticity at the latest time of the simulation, Q(t) ps \u)(X.(ac,t),t)\. 
A vortex line segment L t starting here will intrinsically be "comparable" 
to the maximum of vorticity (as in \u)(K(at,t),t)\ > £7(4)) at late stages 
of the simulation. The assumptions concerning the segment are therefore 
automatically met. In the numerics this procedure is implemented by 
carrying out a precursory identical simulation with a huge number of tracer 
particles (~ 1 million) randomly distributed across the domain. Particles 
that accumulate huge amounts of vorticity are selected for the subsequent 
production run. 

• For the production run, at each instance in time start a vortex line inte- 
gration at X(a,i) along the vorticity direction field. Monitor the max- 
imum curvature ||rt||i«>(i t ) and the maximum vortex line convergence 
11^ • £|U~(L t ) during the integration and calculate X(t). Stop the inte- 
gration, as soon as X(t) reaches a fixed, arbitrary constant C. This defines 
L t . In the numerics this is again implemented with a third-order Runge- 
Kutta integration and bisectioning to obtain the endpoint of L t . 
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Figure 15: Top: Evolution of the length l(t) of the critical vortex filament Lt for different 
Lagrangian fluid elements. The length does not decrease as (T — t) B for any B < 1, which 
would be faster than linear. The Lagrangian collapse of the vortex segment is decelerating 
instead. Bottom: Evolution of the quantities U n and in time. U n does not appear to 
be growing, while U^, though increasing in time, does not exhibit a finite-time blowup as 
l/(T-t). 

• For this vortex line segment L t , calculate the length l(t), and the velocity 
components U n and U^. From the collapse of the length l(t) approximate 
the exponent B. This in turn provides the critical growth exponent A for 
the velocity variables, A cr j t = 1 — B. 

• Compare the increase in U n and £7^ to 1/(T— t) Aclit to distinguish between 
critical and sub-critical growth of velocity. 

This can be interpreted rather intuitively. By prescribing an arbitrarily fixed 
A(i), the vortex line segment is kept relatively geometrically uncritical, as the 
length-scale is always adjusted accordingly. This process of "zooming in" just 
enough to retain the geometric "criticalness" prescribes the rate of collapse to 
a point, at least in the direction of the vortex line. All that is left to check is 
whether the velocity growth in the immediate surrounding is fast enough to be 
compatible with a finite-time singularity. 

The results of the previous section, concerning a point-wise singularity versus 
the blowup of a whole vortex segment, already anticipates that the increase in 
V • £ around the critical vortex line is bounded. If the curvature of the vortex 
line segment remains controllable (which is to be expected from the pictures), 
then just a mild collapse of l(t) occurs. This leaves much room for U n and 
to still be distinguishable from a critical growth. 

Fig. [15] shows the results for the vortex dodecapole initial conditions for a 
fixed constant C . Different choices of the constant C produce identical, rescaled 
results. The top plot pictures the length of the vortex line segment for the 
tracer that is arriving at a position of very huge vorticity at late stages of the 
simulation. The subplot depicts the long-term behavior of the particle entering 
the critical region, while the final stage of length decrease is magnified. The 
decrease in length does not agree with a collapse in final time, but instead the 
shrinkage of the segment decelerates clearly in time. This contradicts a scaling 
in time proportional to (T — t) B for any < B < 1, which would be faster 
than (or, in the limiting case, equal to) linear. It should be noted that for the 
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observed collapse in length, the vortex segment curvature k is the dominating 
term in M(t) = max(||V • $||i~(i t ), ||>t||.E,°°(.E, t )), shadowing the effects of V • £. 
This may lead to a change of regime in the rate of collapse, if V • £ at some 
point exceeds k in quantity. 

It could furthermore be argued that the limit B — > is hard to exclude, 
since the drop in length would be virtually instantaneous in time, with a close 
to constant scaling before. In this limit, the quantities U n and Ut would have 
to grow roughly as 1/(T — t) to still allow formation of a finite-time singularity. 
C/g quantifies the largest difference in axial velocity along the segment. For an 
isolated collapsing vortex tube, this quantity can be expected to not increase 
critically, since the tangential velocity is less likely to rapidly change than the 
radial velocity. However, this initial explanation does not consider the influence 
of the mirror tubes. U„ on the other hand may be interpreted as the velocity 
of the vortex tube's core itself. Again, in an isolated setup this velocity is not 
expected to blow up. 



Fig. 15 (bottom) shows the observed behavior of U„ and E/g in time for the 
tracer as considered above. U n stays roughly constant in time, showing no signs 
of a blowup. E/g, even though increasing in time, does not fit to critical growth, 
in particular not like 1/(T — t) in time. Thus, the assumptions of theorem 2 
are well met. This result therefore poses a strong numerical evidence against a 
finite-time singularity for the class of vortex dodecapole initial conditions. 



5. Conclusion 

In this paper we present numerical evidence against the formation of a finite- 
time singularity for the vortex dodecapole initial condition. We use data ob- 
tained from high resolution adaptively refined numerical simulations to test the 
assumptions presented by geometric blowup criteria. The applied numerical 
method allows for a clearer insight into the formation of the possible singu- 
larity. Most notably, it implies numerical techniques to distinguish between a 
point- wise blowup and the blowup of a whole vortex line segment. Further- 
more, by tracking curvature and spreading of Lagrangian vortex line segments, 
the distinction between singular and non-singular behavior can be made much 
more clearly than the usual approach via BKM . 

In this paper we used vortex dodecapole initial conditions with two different 
vorticity profiles. Comparison of the simulation shows that different vorticity 
profiles yield similar visual and geometrical appearance. This serves as an argu- 
ment that the obtained results may apply to the whole class of vortex dodecapole 
flows. Monitoring the growth rate of Q(t) quantifies the well-known vorticity 
amplification. Amplification by more than two orders of magnitude was reached 
for both vorticity and strain, exceeding by far values achieved by previous sim- 
ulations [40]. Applying this data to BKM would lead to the conclusion that 
a finite-time singularity at time T w 0.72 fits via extrapolation. Yet, as the 
history of Euler simulations has shown, statements obtained by extrapolation 
are to be handled with care. 
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Following the presented argument, a point- wise collapse should coincide with 
a blowup of V • £ at the point of maximum vorticity. As shown, this statement 
can be broadened: A finite-time singularity must either lead to a blowup of V • £ 
at the point of maximum vorticity, or the whole critical vortex line segment 
has to blow up. Utilizing the geometric information obtained via vortex line 
integration from the numerical simulation, it is observed that V • £ does not grow 
in a way to be compatible with a point-wise collapse. Yet, measuring the growth- 
rates on the critical vortex line, high rates of amplification are measured far away 
from the critical region. Even though it is hard to distinguish, whether this 
amplification is critical or sub-critical, this might be interpreted as an evidence 
for the blowup of the complete vortex line segment, even though it suffers exactly 
the same vulnerabilities as extrapolation in BKM. A point-wise blowup, on the 
other hand, cleary contradicts the numerical results up to the time reached. 

Evidence against a blowup of the whole critical vortex line segment is found 
when looking at the geometric properties of Lagrangian vortex line segments. 
It was shown by theorem 2 in [T5] that a blowup of vorticity is directly con- 
nected to the interplay between velocity growth and the collapse of vortex line 
filaments, when maintaining the overall same shape in geometric means (i.e. 
the same X(L t )). Since curvature and V • £ do not increase in order to sup- 
port a finite-time collapse of the segment, velocity components in the vicinity 
of the vortex line filament would have to increase as 1/ (T — t) to support the 
blowup hypothesis. Up to the time reached, this critical velocity growth may 
be excluded by numerical means. This poses a numerical evidence against the 
formation of a singularity in finite time for vortex dodecapole configurations. 
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